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Previously , the performance of the split-symbol moments estimator (SSMP) of 
signal-to-noise ratio (SNR) has been evaluated by means of approximate expressions 
for the estimator mean and variance. These are asymptotic formulas in the sense 
that they become accurate as the number of estimator samples gets large. In the 
present article, exact closed-form expressions are obtained for the same quantities. 
These expressions confirm the accuracy of the previously derived asymptotic results, 
and, unlike the asymptotic formulas , they are useful even when the number of 
samples is small. It is also shown that the conventional split-symbol estimator 
can be trivially scaled to form a signal-to-noise ratio estimator which is precisely 
unbiased (as long as the estimate is based on more than two split-symbols). 


I. Introduction 

The so-called split-symbol moments estimator (SSME) 
of the received symbol signal-to-noise ratio (SNR) was first 
suggested in [1]. In [2], approximate expressions were de- 
rived for the mean and variance of the SSME. These ap- 
proximate expressions were later used as the basis for a 
more detailed performance analysis in [3], which included 
the effects of intersymbol interference as a result of filtering 
the received symbol stream. 

The approximations in [2] result from expanding the 
exact equations for the mean and variance in a series in- 
volving the central moments of the SNR estimator and 
then ignoring all central moments higher than second or- 


der. This type of approximation relies on the fact that 
successively higher central moments of the estimator go to 
zero faster than lower-order central moments as the num- 
ber of samples gets large. 

A natural question is: Under what conditions is it 
valid to throw away the higher-order terms? The results 
quoted for the mean and variance of the SSME (Eqs. 28 
and 31 of [2]) include terms which go to zero as the num- 
ber of sampled split-symbols, n, gets large. These terms 
are proportional to 1/n. How large does n have to be be- 
fore the 1/n terms are guaranteed to dominate the 1/n 2 
terms? Does the answer to this question depend on the 
magnitude of the true SNR? For example, the expression 
for the signal-to-noise ratio of the estimator (Eq. 32 of [2]) 



exhibits peculiar behavior when the true SNR is very small. 
In this case, the equation says that the signal-to-noise ra- 
tio of the estimator actually decreases with n until n gets 
to be about as large as 1/SNR. 

In this article, exact closed-form expressions are ob- 
tained for the mean and variance of of the SSME when the 
symbol stream is unfiltered. These exact results confirm 
the accuracy of the previously known asymptotic expres- 
sions. Unlike the asymptotic formulas, the exact expres- 
sions are useful even when the number of sampled split- 
symbols, n, is small. Furthermore, the basic techniques 
developed for this derivation might also be useful in ana- 
lyzing the more complicated case of filtered data. 


II. Derivation of the Result 

Using the notation of [3], the split-symbol moments 
estimator of SNR is denoted by SNR * and expressed as 


SNR* 
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where m v and m as are average values over n symbols of 
the product and sum squared, respectively, of the outputs 
of half-symbol accumulators, 
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received symbol dj(dj = ±1). The noise samples n fJ are 
assumed to be independent with zero mean and equal vari- 
ance (t 2 /N s . The parameters m and a are the conditional 
mean and standard deviation (given dj\ respectively, of 
the whole-symbol accumulator outputs ( Y a j + Ypj). The 
symbol signal-to-noise ratio SNR is given in terms of m 
and cr by 


(4) 

Conditioned on the sequence of symbol values d = 
(dj, . . . , d n ), the y.j’s are independent Gaussian random 
variables, as are the vectors of half-symbol accumulator 
outputs Y„ = (Yd, . . . ,Y an ) and Yp = (Ypi, . . . , Yp n ). 


The statistics of Y a and Y p are given by 

771 

E {Y aj |d) = E {Ypj |d) = -dj (5a) 

E {Yl |d} = E {Y^ Id} = ^ + y (5b) 

m 2 

E{Y aj Y pj \d} = — (5c) 

m 2 

E {y aj y QJ .|d} = E {YpjYpj, |d} = —djdj. j' # j (5d) 
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E {Y aj Ypy |d} = — dj dj. j ' + j (5e) 
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Vij = y/Sdj + riij 
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i = 0 N, — 1 j = 1,. .. ,n 


Here, N , is the number of samples per symbol and is as- 
sumed to be an even integer, as in [3], and m = N,VS. The 
datum j/,, is assumed to be an unfiltered sample of the jth 


We begin the derivation of our result by rewriting 
Eq. (1) for SNR * in the form 

2 SNR* + 1 = Uj 2 = j- (6) 
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The usefulness of expressing SNR* in the form of Eq. (6) 
results from the fact that U and V are conditionally inde- 
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pendent random variables. Given the data vector d, the 
random vectors u = (iti , . . . , u n ) an d v — ( v i > * • • > v n) are 
jointly Gaussian with 

E {uj |d] = mdj E { u) |d} = m 2 + <r 2 (9a) 
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(for n > 2) 
(for 7i > 4) 
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E{vj\d}=0 E{v]\d} = a 2 

(9b) 

E {uj Vj |d} = 0 = E {uj |d} E {v, |d} 

(9c) 
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E{ujUj' |d} = m 2 djdji = E{u ; |d} E {tij> |d] 
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Equations (9a-e) imply that u and v are uncorrelated 
given d. Because u and v are jointly Gaussian and un- 
correlated, conditioned on d', they are therefore condition- 
ally independent given d. Likewise the random variables 
U and V are conditionally independent given d, because 
they are functions of u and v, respectively [4]. Thus, all 
of the moments of 2 SNR* + 1 can be factored as follows. 

e{(2 SNR * + 1)*} = £d [E {U k \d} E {V' k \d}} (10) 

To evaluate the mean and standard deviation of SNR* 
exactly, it is only necessary to compute the conditional 
moments E{t/|d), E{[/ 2 |d}, E{V -1 |d}. and E{l/" 2 |d}, 
and then calculate the expectation over d. All four of 
these conditional moments can be evaluated in closed form. 
The [/-moments are simply the (positive) moments of a 
noncentral chi-squared random variable. 

E{U\d} = nE{u 7 j \d) =n{m 2 + v 2 ) ( 11 ) 

E (t/ 2 |d) = nE {u 4 |d} + n(n - 1) [E {uj|d}] 2 
= n |(m 4 + &m 2 a 2 + 3<t 4 ) + n(n — 1) (m 2 + a ) 

= 7 i 2 (m 2 -f a 2 ) 2 + n (4m 2 cr 2 + 2a 4 ) (12) 


In Eq. (12), the first equality follows from the indepen- 
dence of the individual uj ’s, and the second equality follows 
from applying the general formula for the fourth moment 
of a nonzero-mean Gaussian random variable [5]. 

The V-moments are the negative moments of a central 
chi-squared random variable with n degrees of freedom. In 
the Appendix, it is shown that 


These V-moments are infinite if n is smaller than the in- 
dicated limits. 


Note that all of the expressions in Eqs. (11), (12), 
(13), and (14) are independent of d. Plugging these re- 
sults’ into Eq. (10) leads to the following expressions for 
the moments of the estimator SNR (assuming n > 2 or 
7 i > 4 as appropriate): 


E{2SNR* + l} = 


n (m 2 + a 2 ) 
(n - 2 )<t 2 


= (2 SNR + 1)^2 
(15) 


E |(2 SNR* + l) 2 } = (m 2 + cr 2 ) 2 + 1 (4m 2 <r 2 + 2<r 4 ) 


(n - 2)(n - 4)<r 4 

= J^(2 SNR + l) 2 + i (8 SNR + 2) 

x (16) 

(n - 2)(n - 4) 

The mean and variance of SNR* are calculated from 
Eqs. (15) and (16) as 

E {SNR* } = i [E {2 SNR* + 1} - 1] 

= SNR+ — 5-;r(2 SNR + 1) (for n >2) 
n — 2 

(17a) 

var {SNR*} = ^var {2 SNR* + 1} 

= I [jE {(2 SNR* + l) 2 } - (E {2 SNR* -1- l}) 2 ] 


(for n > 4) 


(17b) 
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In comparison, in [2] the following approximate asymptotic 
results were obtained: 


E {SNR?* } = SNR (for n > 2) (20a) 


E{SNR*} « SNR + i(2 SNR+ 1) (I8a) 

var{SAT iT} « I [2 SNR 2 + 4 SNR+ l] (18b) 

By comparing Eq. (17) with Eq. (18) it is seen that the 
asymptotic expressions of [2] are accurate if , and 

can all be approximated by 1. 

In addition to the mean and variance, higher-order 
moments of the SSME may also be evaluated by applying 
known formulas for the moments of noncentral chi-squared 
random variables ({/-moments) and for the moments of 
the reciprocal of a central chi-squared random variable 
(V-moments), and then solving Eq. (10) algebraically for 
the corresponding moments of SNR* . 

Finally, it is seen from Eq. (15) that it is easy to con- 
vert SNR* into an unbiased estimator SNR** (for n > 2) 
by making the definition 

2 SNR** + 1 = (2 SNR* + 1) (19a) 

i.e., 

SNR** = SNR * - - (2 SNR* + 1) (19b) 

n 

Then, assuming n > 2 for Eq. (20a) and n > 4 for 
Eq. (20b), 


var {SNR**} = ^^2— 2 SNR 2 + (4 SNR + 1) 

(for n > 4) (20b) 

This should be the preferred form for the split-symbol mo- 
ments estimator because, in addition to being unbiased, 
SNR** also has a slightly smaller variance than SNR* . 


III. Conclusion 

The technique introduced in this article enables one 
to compute exact closed-form expressions for the perfor- 
mance of the split-symbol moments estimator, for the case 
of unfiltered, undegraded data. The essential trick required 
for this computation was to algebraically manipulate the 
estimator formula into the form of a quotient of condi- 
tionally independent random variables. The calculations 
in this article confirm the accuracy of previously derived 
asymptotic expressions for unfiltered data. In addition, the 
exact performance expressions, unlike the asymptotic for- 
mulas, are useful even when the number of sampled split- 
symbols, 7 i, is small. The exact formulas show that as 
long as n > 2 the conventional split-symbol estimator can 
be trivially scaled to form a signal-to-noise ratio estima- 
tor which is precisely unbiased. Furthermore, the same 
techniques developed here may be applied to obtain sim- 
plified .expressions (though not closed form) for the more 
complicated case when the data is filtered. 
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Appendix 

The Moments of the Reciprocal of a Chi-Squared Random Variable 


Conditioned on d, the random variable V is a cen- 
tral chi-squared random variable with n degrees of free- 
dom. According to [6], its conditional probability function 
is therefore given by 


p(K|d) = 


v {n/2)-l e -V/2v* 

(2<r 2 )" /2 r(n/2) 


V > 0 


(A-l) 


too y(n/ 2 )-fc-l e -V 72 < 7 3 

^ / d\f 

Jo (2^) {n/2) - k r[(n/2) - k] 

(r> 2\~ k r[( n /2) - k] jf ^ 0/. 

- ^ } T(n/2) 11 n> 


(n - 2)(n - 4) • • • (n - 2 k)a 


2k 


(A-2) 


where T(-) is the gamma function. Thus, the conditional 
(negative) moments of V can be calculated (for k > 0) as 


roo y(nl2)-k-\ e -VI2<, 2 

^ ^ “ Jo (2cr 2 ) n/2 iYn 


(2cr 2 ) n/ ^ T(n/2) 


dV 


= M 


-k r[(n/2) - fc] 
T(n/2) 


The next-to-last equality in Eq. (A-2) follows from noting 
that, if n > 2k, the rescaled integral in the preceding line 
is simply the integral of the probability density function 
of a chi-squared random variable with n — 2k degrees of 
freedom. The last equality in Eq. (A-2) follows from the 
gamma function recurrence formula, T(z + 1) = zY(z). 
Finally, it is easy to see that the integral in Eq. (A-2), 
and hence the corresponding moment of V , is infinite if 
n < 2k, because in this case yC"/ 2 )-*- 1 becomes infinite 
near V = 0 at least as fast as V~ l . 
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